function[Xind,Yind,ie]=voronoi_drawxy(width,height,equation_set,equationsymbol_set,relation_slack)
        [X,Y]=meshgrid(10:5:width,10:5:height);%10:20
        ind=(1==1);
        for i=1:size(equation_set,1)
            if relation_slack(i,1)==1
                if equationsymbol_set(i,1)==1
                    ind=ind&(equation_set(i,1)*X+equation_set(i,2)*Y+equation_set(i,3)>=0);
                end
                if equationsymbol_set(i,1)==-1
                    ind=ind&(equation_set(i,1)*X+equation_set(i,2)*Y+equation_set(i,3)<=0);
                end
            end
        end
        ie=isempty(X(ind));
        Xind=X(ind);
        Yind=Y(ind);
end